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i"ff1etodo per La determinazione delta superficie tridimensionale di un oggetto" | 
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L BIASSOHTO 



La presente invenzione si riferisce ad un metodo per la detenninazione della superficie tridimensionale di un 
oggetto. 

In una sua forma di realizzazione il metodo per la determinazione della superficie tridimensionale di xm oggetto 
comprende le fasi di: definire (1) le coordinate di una pluralita di punti di detto oggetto; definite (2) xma matrice 
tridimensionale di celle che contenga detto oggetto alle quali e possibile associare un vSore; determinare (3) la distanza 
tra ogni centro di dette celle di detta matrice tridimensionale di celle ed il pimto piu vicino di detta pluralita di punti di 
detto oggetto; fissare (4) il valore di alcune celle di detta matrice tridimensionale di celle ad un primo valore prefissato; 
determinare (7) il valore che assume ogni cella di detta matrice tridimensionale di celle, ad eccezione di dette alcune 
celle, mediante la seguente formula 

m.^ +1) = -T—^ ! 

dove^ Xi rappresenta le coordinate del centro della cella i_esima, 

F{x^ , f) rappresenta il valore della cella i_esima al tempo ; 

rappresenta detta distanza, j 

"W rappresenta un secondo valore prefissato, e j 

j indica un intomo di celle della cella i_esima; j 
determinare (9) la somma in modulo delle variazioni del valore di ogni cella tra il tempo / ed il tempo ^+7; •' 

ripetere (10) detta fase di determinare il valore che assiune ogni cella di detta matrice tridimensionale di celle se detta } 

somma h maggioie di un terzo valore prefissato. (Fig. 1). 
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DESCRIZIONE 
dell'invenzione indvistriale avente per titolo: 

"Metodo per la determinazione della superficie tridimensionale di un 
oggetto". 

a riome: POLITECNICO DI MDLANO 



Hp 4* 



La presente invenzione si riferisce ad un metodo per la determinazione 
della superficie tridimensionale di un oggetto e relative progranama per 
calcolatore. 

Per descrivere ed analizzare la geometria di oggetti che devono essere 
realizzati sono noti metodi di modellazione che utilizzano calcolatori. In 
particolare, quando si vuole ottenere una superficie tridimensionale chiusa da 
un insieme di punti, acquisiti da im oggetto reale, si fa normalmente uso di 
una tecnica denominata level-set che descrive una superficie chiusa come una 
superficie di livello di una funzione volumetrica (o funzione implicita). 

La tecnica level-set prevede Tevoluzione temporale della funzione 
volumetrica, secondo una opportuna equazione alle derivate parziali, 
tipicamente appartenente alia famiglia di equazioni di Handlton - Jacobi, 

La funzione volumetrica descrive la distanza di ciascun punto dalla 
superficie chiusa e tale distanza sara dotata di segno positivo o negativo a 
seconda che il punto considerato sia intemo o estemo alia superficie. Quindi 
rinsieme di punti che possiedono una distanza nulla dalla superficie S 
rappresentano la superficie stessa. La tecnica prevede di far evolvere la 
superficie chiusa, definita dall'equazione di Hamilton - Jacobi in modo che 
essa circondi ed aderisca alia nuvola di punti 3D che descrive I'oggetto. 
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Viene utilizzata qviesta tecnica in quanto in grado di adeguarsi alle varie 
tipologie di oggetti ed alle pit. disparate tecniche di acquisizione dei punti. 
Tale tecnica infatti, non richiede alcuna informazione aggiuntiva sulla 
topologia dei punti acquisiti ed h insensibile alia tecnica utilizzata: la densitsi 
dei punti pu6 variare nelle diverse zone dell'oggetto considerate; non h 
richiesta alcuna informazione sulla sequenza di acquisizione dei punti o sulle 
loro relative posizioni spaziali; ed il numero di oggetti chiusi separati tra loro 
b ininfluente per Talgoritmo. 

Poich6 questa tecnica richiede Taggiomamento di una funzione 
volumetrica ad ogni passo di evoluzione del fronte, il costo computazionale e 
solitamente molto elevato, al pxanto da rendere il metodo di scarsa utility 
pratica. Sono noti metodi che consentono di ridurre il costo computazionale di 
questa tecnica, ad esempio limitandosi ad aggiomare la funzione volumetrica 
in una zona circoscritta circostante il fronte in evoluzione. cid nonostante il 
costo computazionale rimane proibitivo. 

Sono anche note tecniche alternative a quelle del level-set Ad esempio, 
sono noti metodi come quelli della triangolazione di Delaunay, del NURBS 
(Non Uniform Rational B-Splines) e del HREF (Hierarchical Radial Basis 
Function). Ciascimo di questi metodi ha delle limitazioni ed impongono una 
serie di vincoli sull'insieme dei punti acquisiti, come ad esempio hanno una 
forte dipendenza della complessita computazionale dal numero di punti, 
oppure hanno la necessita di avere una distribuzione spaaale dei punti 
acquisiti quanto piii uniforme possibile. 

In vista dello stato della tecnica descritto, scopo della presente invenzione 
e quello di prowedere ad un metodo per la determinazione della superficie 
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tridimexisionale di im oggetto veloce e che non abbia gli inconvegnenti della 
arte nota. 

In un suo primo aspetto la presente invenzione si riferisce ad un metodo 
per la determinazione della superficie tridimensionale di un- oggetto 
coinprendente le fasi di: definire le coordinate di una pluralita di punti di detto 
oggetto; definire una matrice tridimensionale di celle che contenga detto 
oggetto alle quali 6 possibile associare un valore; determinare la distanza tra 
ogni centre di dette celle di detta matrice tridimensionale di celle ed il punto 
piu vicino di detta plurality di punti di detto oggetto; fissare il valore di alcune 
celle di detta matrice tridimensionale di celle ad un primo valore prefissato; 
determinare il valore che assume ogni cella di detta matrice tridimensionale di 
celle, ad eccezione di dette alcime celle, mediante la seguente formula 

F{x,J^ 1) = ^ 

J 

dove 

x,rappresenta le coordinate del centro della cella i_esima, 
^(^/»0 rappresenta il valore della cella i_esima al tempo 
. V, rappresenta detta distanza, 
w rappresenta un secondo valore prefissato, e 
j indica un intomo di celle della cella i_esima; 

determinare la sonuna in modulo delle variazioni del valore di ogni 
cella tra il tempo / ed il tempo r+i; ripetere detta fase di determinare il valore 
che assume ogni cella di detta matrice tridimensionale di celle se detta somma 
h maggiore di un terzo valore prefissato. 
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In un suo secondo aspetto la presente invenzione si riferisce ad un 
programma per calcolatore comprendente un codice di programma che esegue 
tutte le fasi del metodo per la determinazione della superficie tridimensionale 
di un oggetto quando detto programma e eseguito su detto calcolatore. 

In un suo terzo aspetto la presente invenzione si riferisce ad xin 
prQgramma per calcolatore registrato su un supporto utilizzabile da detto 
calcolatore per controUare Tesecuzione di tutte le fasi del metodo per la 
determinazione della superficie tridimensionale di vm oggetto. 

Grazie alia presente invenzione, sostituendo Tequazione differenziale alle 
derivate parziali nomialmente utilizzata in letteratura (Hamilton- Jacobi) con 
ima che si presta molto bene a descrivere funzioni discontinue quali sono 
quelle che descrivono Tinterfaccia fra due fluidi immiscibili, ossia una delle 
equazioni di Navier-Stokes, derivata dalla fluidodinamica, si riducel 
notevolmente il costo computazionale, ed 6 possibile ottenere una 
convergenza molto pitL rapida. La superficie pud essere inteipretata come il 
fronte tra due fluidi viscosi che si muove fino ad awiluppare completamente 
I'oggetto. In particolare, agendo sul valore di viscositi dei fluidi impiegati 
risulta inoltre possibile controUare il livello di rugosity consentendo di 
ottenere superfici Usee anche qualora i punti presentassero deUe deviazioni 
dal loro corretto posizionamento spaziale. Inoltre, il metodo proposto 
consente di introdurre elementi e strategic che consentono un controUo fine 
deU'evoluzione del fronte, basato su principi di sempUce interpretabilit^ 
fisica. Ad esempio, e ppssibUe introdurre una turbolenza artificiale nei fluidi 
che agevoU la penetrazione aU'intemo di cavita particolannente strette. 
Risulta cosi possibUe fomire una determinazione corretta anche di 
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configurazioni di punti che tradizionalmente erano state negate. L'equazione 
di Navier-Stokes, grazie alia sua flessibiliti, puo essere impiegata per 
determinare piu funzioni volumetriche che descrivono diverse parti 
dell'oggetto considerate. Oggetti particolannente complessi o con elevate 
risoluzioni possono qnindi essere rappresentati calcolando separatamente 
revoluzione della funzione volumetrica snlle diverse parti che li 
compongono. I singoli risultati ottenuti vengono poi ricomposti per ottenere la 
funzione volumetrica complessiva impiegando un semplice operatore 
matematico. 

Le caratteristiche ed i vantaggi della presente invenzione risulteranno 
evidenti dalla seguente descrizione dettagliata di una sua forma di 
realizzazione pratica, illustrata a titolo di esempio non limitativo neirunito 
disegno, nel quale la figura 1 mostra un diagramma di flusso del metodo per 
la determinazione della superficie tridimensionale di un oggetto. 

La presente invenzione si basa sulFutilizzo della equazione dififerenziale 
alle derivate parziali di Navier — Stokes per la conservazione della massa nella 
quale viene utilizzata una ridefinizione del vettore velocita secondo un criterio 
descritto nel seguito. Essa descrive il moto di due fluidi immiscibili che 
convergono, I'xmo dairintemo e Taltro dall'estemo, verso la nuvola di punti 
che descrive Toggetto considerate. Al termine deU'evoluzione, quando 
entrambi i fluidi aderiscono ai punti dell*oggetto, la superficie di separeizione 
tra i due coincide con la superficie dell'oggetto stesso. 

La funzione volumetrica descrive il contenuto dei due fluidi in ciascuna 
cella elementare o voxel, dove con voxel si intende Telemento elementare 
della griglia volumetrica, cosi come il pixel Adene utilizzato per una griglia 



Dr. Ing. Enrico Mittler 



bidimensionale. Per distinguere i due fluidi abbiamo attribuito ad essi massa 
di segno opposto. Pertanto una cella nella quale la funzione volumetrica valga 
-1 contiene esclusivamente il fluido intemo, mentre un valore di +1 indica che 
e presente esclusivamente il fluido estemo. 

- Se tra due celle attigue la funzione volumetrica cambia di segno implica 
che tra essi passa la superficie di separazione e quindi anche la superiHcie 
deiroggetto- L'esatto punto di passaggio viene determinato mediante 
operazioni di interpolazione. 

La posizione esatta della superficie di contatto tra i due fluidi viene 
quindi identificata dalla posizione del level-set zero (ossia Tinsieme di punti 
con distanza nulla dalla superficie). 

La legge per la conservazione della massa 6 uno dei principi fondamentali 
della fisica classica ed e indipendente dalla natura dei fluidi e dalle forze che 
agiscono su di essi; essa esprime il principio secondo il quale in un sistema 
fluidodinamico la massa non pud scomparire o venir creata se non in presenza 
di pozzi o sorgenti. 

II valore F della funzione volumetrica in ciascuna cella indica la quantity 
di fluido al suo intemo e tale valore viene modificato dalle forze agenti sul 
fluido, dalla sua velocita, dalla diffiisione proveniente dalle celle attigue e 
dalle sorgenti o pozzi. 

L'evoluzione di un fluido all'intemo di un sistema puo essere descritta 
tramite il vettore di flusso G che esprime la quantita di fluido che attraversa 
una superficie in un dato intervallo di tempo. 

Nella sua formulazione classica la legge di conservazione per la massa 
afierma che la variazione per xmita di tempo della quantity di materia F 
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airintemo di un volume /2, 




e pari al contributo netto dei flusso che attraversa la superficie estema S 
di al contributo Qy di sorgenti (o pozzi) di volume al suo intemo, ed al 
contributo di sorgenti 2^ orientate, lungo la sua superficie S. I contributi 
dovuti alle diverse sorgenti possono quindi essere scritti come: 



La fomiulazione generate per Tequazione di conservazione della massa 
puo quindi essere scritta come: 



II segno negative davanti alFintegrale di flusso e dovuto al fatto che 
Fareola dSb orientata verso Testemo di /2 

Utilizzando il teorema di Gauss e possibile riscrivere la precedente 
equazione utilizzando esclusivamente termini di volume: 



II vettore densita di flusso G puo essere scomposto in uria componente 
diffiisiva G^, ed in una convettiva G^ , la parte convettiva del vettore di flusso 
descrive il fenomeno del trasporto dovuto alle forze esteme e viene definita 
come prodotto tra la velocita v del flusso e la quantita di materia trasportata F: 
risulta pertanto Gc-vF, 

La componente diffiisiva e invece derivante dal rimescolamento, presente 
anche nei fluidi a riposo, dovuto aU'agitazione termica, ed h solitamente 
proporzionale, secondo la legge di Pick, al gradiente di F 



n s 
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^ys/F dovG 6 la costante di diffusivita. L'evoluzione temporale 
della funzione volximetxica F dovuta alia difEusiviti h regolata (come risulta 
verificabile imponendo la conservazione della massa) dalFequazione 
differenziale alle derivate parziali del calore e quindi il suo contributo risulta - f^jj^^T; 

principalmente connesso con una regolarizzazione ed un appiattimento della 
funzione F. L'effetto del contributo diffusivo h quindi stato preso in lb 
considerazione applicando un filtraggio gaussiano alPintera funzione 
volumetrica. 

L'intensita e I'orientazione della velocity v del fluido e determiriata in 
ciascuna cella della matrice. Tale velocita e stata definita orientando il flusso 
verso il punto della nuvola 3D (che descrive I'oggetto) piu vicina alia cella, ed 
il modulo di tale vettore e stato scelto direttamente proporzionale alia 
distanza: in tal modo entrambi i fluidi rallentano progressivamente 
awicinandosi alia nuvola di punti e convergono alia superficie di conta 
finale. L'equazione che e stata quindi utilizzata per descrivere Tevoluzio 
dei due fluidi b la seguente: 





dF(Xi,t) _ ^ 

at 



Ji5'(5c,0|v(jc)|-e ^ dx 



J|v(jc)|.e ^ 



La variazione locale della quantita di fluido contenuto nella cella ijesima 
F(x,,Otra due istanti successivi risulta pertanto determinata dal contributo 
convettivo derivante da ciascuna cella F(x,t)\v(x)\ pesato con una gaussiana 
la cui deviazione standard o descrive il grado di diflfiisione del fluido tra le 
diverse celle. 



-9- 
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La velocita v(x)h stata definita per ciascuna cella proporaonale alia 
distanza tra U centre della cella ed il piti vicino punto appartenente alia nuvola 
3D. In particolare il modulo di tale velocita e pari a |x, -^"dove x, h H 
centro della cella i_esima, mentre p h'H punto ad esso pit vicino della nuvola 
di punti (il punto a distanza minore dal centro della cella i_esima), a e invece 
un parametro che regola I'andamento della funzione velocity 
Sperimentalmente e stato detenninato che valori accettabili per a sono 
compresi tra 1,5 e 2,5. Tanto maggiore h tale valore tanto pit la velocita sar^ 
lenta in prossimiti dei punti mentre sack rapida a grande distanza. Tale 
caratteristica, anche se molto importante per una evoluzione rapida dei valori 
della matrice pu6 comunque portare ad una convergenza finale lenta. Un 
valore preferito del valore di a 6 pari a 2.0. 

Un'altra propriet^l dei fluidi utilizzata nel modello 6 la viscosity il cui 
valore determina la deviazione standard <x all'intemo dell'equazione. Tale 
propriety descrive I'attrito viscoso che agisce intemamente al flxiido 
limitandone la mobilitSi. Considerando tale propriety impediamo che i fluidi 
attraversino la superficie ricostruita fluendo tra gli interstizi dei diversi punti. 
Questa propriety 6 particolarmente importante in quanto tramite essa 6 
possibile onoraie correttamente particolari disposizioni di punti quali lame 
sottili od oggetti particolarmente appimtiti. Grazie alia viscosita, i fluidi non 
attraversano la superficie la quale resta aderente alia lama o alia punta. Risulta 
inoltre possibile ridurre la curvatura dell'intera superficie ottenendo supeifici 
particolarmente levigate e riducendo la rugosita superficiale che pu6 anche 
essere origmata dall'errore di misura. Occorre poro ricordare che valori 
elevati possono comunque limitare la capacity dell'algoritmo di rappresentare 
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oggetti con superfici particolannente frastagliate e di consentire Tafflusso del 
liquido all'intemo di zone concave. Per quanta sopra detto la viscosity h un 
parametro la cm scelta viene determinata a seconda della tipologia 
deU'oggetto considerate, dall'imprecisione del sistema di acquisizione 
utilizzato e dal grado di fedelt^l col quale desidera che la superficie finale 
aderisca alia nuvola di punti. 

Un ulteriore aspetto che e state utilizzate dalla modellizzazione fisica dei 
fluidi 6 la turbolenza, la quale ixatta la presenza di gorghi che descrivono 
variazipni locali della ftmzione velocity. 

La presenza della turbolenza pu6 sviluppare localmente pressioni elevate 
e possono generare instability nel flusso ma alio stesso tempo puo favorire la 
penetrazione del fluido all'intemo di aperture molto piccole ed impedime il 
ristagno in zone con velocita modeste. 

La presenza di vortici e stata introdotta modulando il campo di velocity 
con piccole oscillazioni nelle diverse direzioni attomo alia configurazione 
iniziale. A causa della stretta correlazione tra il vettore velocita in ciascuna 
cella ed il vettore che punta verso il punto piu vicino della nuvola 3D, le 
oscillazioni locali del campo di velocita possono anche essere interpretate 
come una deformazione elastica della nuvola di punti nelle varie direzioni la 
quale puo agevolare I'ingresso del fluido all'intemo di zone anguste. 

La modulazione della fimzione velocity e applicata ciclicamente lungo i 
tre assi spaziali moMplicando per un fattore i termini v delle velocita delle 
celle che giacciono orientati parallelamente a tali assi rispetto alia cella 
i_esima in esame e la sua intensita e controllabile estemamente, in pardcolare 
ad ogni iterazione viene amplificata la velocita in una data direzione e/o 
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ridotta di una pari quantita la velocM nelle altre due direzioni. Owianaente 
maggiori saranno i coefficienti di espansione utilizzati maggiore sari la 
turbolenza. Valori troppo elevati possono per6 determinare un eccessivo 
rimescolamento anche tra i due fluidi suUa superficie di contatto inducendo la 
fuoriuscita del fluido intemo o un'eccessiva penetrazione del fluido estemo 
neiroggetto considerato. 

L'equazione precedente che h stata utilizzata per descrivere Tevoluzione 
dei due fluidi h stata discretizzata per poter eseguire le operazioni piu 
facilmente mediante un calcolatore ed e descritta dalla seguente equazione. 

F(x„r + 1)-F(x„0 = v^^w^v, ""^"-'^ 

J 

dove 

jc^rappresenta le coordinate spaziali (array tridimensionale) della cella 
iN^^esima, 

F(Xf,t) rappresenta il valore della cella i^esima al tempo 
V,. rappresenta la velocita (o la distanza come precedentemente 
definita), 

w rappresenta un valore prefissato di viscosita, e 
j 6 un valore prefissato che indica im intomo di celle della cella 
i_esima. 

In particolare, la sommatoria indice j considera il contributo delle celle 
immediatamente piii vicine al punto in esame. Preferibilmente, e stato scelto 
di limitare tale interazione esclusivamente alle prime 6 celle vicine a ciascuna 
cella in esame, ossia che condividono con esso un'intera &ccia della cella 
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cubica (le celle a distanza 1 dalla cella in esame). 

Ci riferiamo ora alia figura 1 che mostra il diagramma di flusso del 
metodo per la determinazione della superiBcie tridimensionale di un oggetto. 
Tale diagramma di flusso descrive come un programma per calcolatore 
effettua il calcolo per la deteraiinazione della superficie tridimensionale di un 
oggetto. 

Fase 1 - Definizione dati. X2^^^^><v 
In questa fase vengono definite le coordinate di una pluralit^L 
deiroggetto che si vuole rappresentare. L'insieme dei pimti viene des^tto 
come un insieme di vettori a 3 componenti, contenenti le coordinate spa 
di ciascun punto. Non viene richiesta nessuna informazione aggixmtiva 
relativa ai punti. 

Fase 2 - Definizione matrice di celle. 

In questa fase viene descritta la matrice di celle (o griglia>^j 

tridimensionale) che definisce il dominio di evoluzione del valore di F. nT^ 

V % 

numero di celle 6 scelto in base aUe dimensioni dell'oggetto che si vuoie. . 
rappresentare. Si determina il numero di celle che verranno utilizzate in base\<{/7 
airestensioni lungo i tre assi deU'oggetto. Mantenendo quindi una forma 
cubica per ogni cella elementare della matrice, viene definita una griglia che 
contiene esattamente il nostro oggetto. Risulta anche possibile definire una 
forma delle celle non cubica, ossia si puo impostare una diversa risoluzione 
della matrice lungo i diversi assi. Questo puo essere molto importante per 
oggetti nei quali sia necessaria una elevata risoluzione soltanto in una 
dimensione come ad esempio bassorilievi, facciate di edifici, ecc. 

E' anche possibile definire una matrice ruotata rispetto £^li assi secondo 
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cui b stato acquisito Toggetto ed, in particolare, allinearla con le direzioni 
principali deiroggetto stesso. Risulta cosi possibile minimizzare il numero di 
celle impiegate pur mantenendo la medesima risoluzione. Ad ogni cella 
appartiene irn centro con le appropriate coordinate, ed ad ogni cella sar^ 
possibile associare un valore 

Fase 3 - Determinazione della velocity 

In questa fase viene definito il campo delle velocity che competono a 
ciascuna cella della griglisu Come abbiamo visto precedentemente il campo di 
velocitii h strettamente dipendente dalla distanza di ciascun centro delle celle 
dal punto piu vicino della nuvola di punti. Quindi, il campo delle distanze e 
una matrice tridimensionale delle medesime dimensioni della griglia. Per 
ciascuna cella viene determinato il pxmto della nuvola 3D che gli risulta piu 
vicino (a distanza minore) e viene calcolata la distanza come \xf -p^ tra tale 
punto ed il centro della cella. Tale valore viene memorizzato nella cella 
equivalente della matrice delle distanze. 

Si noti che non e necessario conoscere la direzione delle velocita ma 
soltanto il modulo. 

E' importante notare come, la dipendenza della complessita 
compxitazionale di questo metodo dal numero di punti della nuvola 3D sia 
limitata esclusivamente alia determinazione del campo di velocita, una sola 
volta airinizio della procedura, e tutta la fase evolutiva 6 completamente 
indipendente dal numero di punti campione utilizzati. 

Fase 4 - Disposizione delle sorgenti. 

Sono possibili diverse disposizioni delle sorgenti. Una prima disposizione 
prevede che le sorgenti esteme siano poste esclusivamente sui bordi della 
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matrice (i quali verranno quindi sempre mantenuti al valore +1 durante tutta 
I'evoluzione) e tutti le restanti ceUe siano invece riempiti di fluido estemo 
(valore pari a -1). Una seconda disposizione invece prevede che vengano 
fomiti anche una o piii sorgenti interne, in tal caso le celle del contomo 
saranno sempre mantenuti a +1, le celle interne avranno valore pari a 0, 
mentre le sorgenti interne saranno sempre mantenute a -1. Questo metodo 
consente un'evoluzione ancora piu veloce ma richiede delle informazioni 
iniziali relative alia collocazione delle sorgenti interne j^iuntive che sono 
dif&cilmente disponibili. 

Fase 5 - Fissaggio parametri di viscosity e turbolenza. 
La viscosit^i (nella formula discretizzata rappresentata da w) regola la 
diffusione del fluido e quindi il peso dell'interazione tra ogni singola cella ed i 
suoi vicini. II valore preferito di w risulta compreso tra 0,1 e 0,9. 

I parametri relativi alia turbolenza indicano come le velocity del fluido, in 
ogni cella, vengano alterate ad ogni iterazione al fine di simulare I'effetto 
turbolento. In particolare risulta possibile fissare per ciascimo dei 3 assi il 
coefl&ciente per il quale weak moltiplicata ciclicamente la velocity in tale 
direzione. I fattori moltiplicativi per la velocita variano preferibilmente da un 
minimo di 0,5 fino ad un massimo di 1,5. L'uso del regime turbolento h 
facoltadvo. 

Fase 6 - Awio della procedura. 

Per velocizzare le operazioni, ed in modo facoltativo, in questa fase viene 
predisposta ima ulteriore matrice tridimensionale delle dimensioni della 
matrice delle celle nella quale viene contenuta I'informazione 
sull'aggiomabilita o meno delle diverse celle. Se una cella t aggiomabile il 
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suo contenuto sar^ processato dairalgoritmo evolutive e quindi il contenuto 
di liquido al suo intemo potr^ variare tra due iterazioni. Celle non aggiomabili 
sono le sorgenti (che vengono sempre mantenute al medesimo valore, +1 o - 
1). Nel corso deli'evoluzione possono diventare non aggiomabili anche quelle 
celle per i quali la variazione di liquido tra due iterazioni non e stata 
significativa (al di sotto di una certa soglia impostata) e qiiindi tali celle 
vengono congelate fintanto che il fronte evolutive non giunga a lambirli. 
Fase 7 - Aggiomamento celle. 

E' la fase principale del metodo nella qxiale si determina il valore della 
quantita di fluido presente in ogni cella della matrice tridimensionale di celle. 

L'impiego della gaussiana airintemo dell'integrgle prevedrebbe FanaUsi 
ad ogni iterazione deirinterazione di ciascuna cella con tutto la matrice delle 
celle, il che risulterebbe particolarmente gravoso in termini computazionali. 

Nel caso con la formula discretizzata, risulta possibile e preferibile 
limitare Testensione di tale interazione esclusivamente ad un intomo di celle 
piix vicine alia cella i^esima, che possono essere ad esempio 6 o 18 o 26 (la 
cella considerata e al centro di una matrice 3x3x3), cio6 quelle che si trovano 
rispettivamente a contatto di faccia, di spigolo o vertice (ossia a distanza 1, ^2 
e >/3 dalla cella in esame, nel caso di celle cubiche) daUa cella i^esima. 
Preferibilmente sono state scelte le 6 celle piii vicine, ossia le 6 celle che 
contattano le 6 facce della cella i^esima, 

Una volta aggiomato la cella, viene valutata la variaaone rispetto al 
valore precedente e se tale valore e mferiore ad una soglia impostata 
dall'utente il contenuto della cella viene congelato ed la cella viene dichiarato 
non aggiomabile. 
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Fase 8 - Vsualizzazione 

Risulta possibile, se di interesse, visualizzare il level-set zero, ossia la 

superficie di separazione tra i due fluidi. L'algoritmo utilizzato e quello 

comunemente noto col nome di Marching Cubes (si veda ad esempio 

rarticolo "Marching cubes: a high resolution 3D surface construction 

algorithm", di William E. Lorensen e Harvey E. Cline, da Computer Graphics, 

Volume 21, Nxmiero 4, Luglio 1987) mediante il quale 6 possibile ottenere 

un'efficace triangolazione del level-set zero ed una risoluzione sub-pixel/^ 

grazie all'interpolazione lineare effettuata su celle contigue. La triangolazilne*;^ 

cosi ottenuta viene rappresentata a video, ad esempio, tranute un programma 

di interfaccia denominato OpenGL che consente una visualizzazione in tempo f^^^^^J^^ 

reale ed una interazione diretta con Tutente il quale pu6 interatdvamente 

cambiare a suo piacimento il punto di vista analizzando Tevoluzione 

l^j o : I 10,33 lEftrb^ 

level-set in tempo reale. ^ 
Fase 9 - Determinazione delle variazioni di flusso. \ s 

Al terinine di ogni evoluzione viene valutata la quantita complessiva di 
fluido che h stata mossa. 

Fase 10 - Evoluzione ciclo. 

Se la quantity complessiva di fluido che e stata mossa e inferiore di una 
soglia prefissata Tevoluzione viene interrotta alia fase 11, altrimenti effettua 
un'ulteriore iterazione. 

Alia fese 11, quando si e tenninato il calcolo, si ottiene una matrice 
tridimensionale di punti che definisce i valori che assume ogni punto la 
fiinzione volumetrica. 

Risulta possibile applicare un filtro passa basso a tale matrice in modo da 
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ottenere una superficie ricostruita piu levigata rimuovendo possibili errori di 
campionamento senza agire sui parametri di viscosita (che potrebbero 
pregiudicare la corretta determinazione di zone concave). L'implementazione 
di tale filtro b analoga a quella per effettuare un filtraggio mono e 
bidimensionale: I'intera funzione volumetrica viene convoluta con una matrice 
tridimensionale contenente i valori del filtxo. Ad esempio h possibile 
realizzare in filtro passa basso impiegando una matrice 3x3x3 nella quale 
Telemento centrale ha il valore massimo e gli altri elementi possiedono un 
valore dato dal valore di una funzione gaussiana dipendente dalla loro 
distanza del centro della matrice. 

Con il presente metodo, nel caso in cui risulta possibile e conveniente 
scomporre una superficie di grosse dimensioni in parti distinte, risulta 
possibile calcolare Tevoluzione singolarmente su ciascuna parte e ricomporre 
poi le varie funzioni volumetriche. In questo modo risulta anche possibile 
utilizzare risoluzioni differenti per le diverse parti per poi ricomporle su 
un'unica matrice comune. E' preferibile, al fine di migliorare il raccordo tra le 
diverse parti, che tra le diverse zone definite^ vi sia sovrapposizione, ossia che 
parte dei pimti estremi appartenenti ad una zona vengano inglobati anche nella 
zona contigua. Una volta che si hanno le matrici delle varie zone, esse 
vengono disposte su ima matrice comime con risoluzione pari alia massima 
risoluzione impiegata. Nel caso in cui un il contenuto di ima cella della griglia 
finale venga contemporaneamente descritto da due fiuizioni volxmietriche (in 
quanto si sovrappongono) viene scelto come valore per tale cella il minimo tra 
tutte le fiinzioni volumetriche che lo descrivono. 

II metodo per la determinazione della superficie tridimensionale di un 
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oggetto, qui descritto, h trascrivibile in un codice di progranmia, in modo ben 
noto da un tecnico del ramo, che pu6 essere memorizzato su qualsiasi tipo di 
memoria o supporto (floppy, CD) e/o reso eseguibile da un calcolatore. 
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RIVENDICAZIONI 
1 . Metodo per la deteniainazione della superficie tridimensionale di un 
oggetto comprendente le fasi di: 

definire (1) le coordinate di una pluralita di punti di detto oggetto; 
' definire (2) una matrice tridimensionale di celle che contenga detto 
oggetto alle quali h possibile associare un valore; 

determinare (3) la distanza tra ogni centre di dette celle di detta matrice 
tridimensionale di celle ed il punto piu vicino di detta pluralita di punti di 
detto oggetto; 

fissare (4) il valore di alcune celle di detta matrice tridimensionale di 
celle ad un primo valore prefissato; 

determinare (7) il valore che assmne ogni cella di detta matrice 
tridimensionale di celle, ad eccezione di dette alcxme celle, mediante la 
seguente formula 

F(x,,t^D= 

J 

dove 

3C|rappresenta le coordinate del centre della cella ijesima, 
F(3c; , 0 rappresenta il valore della cella i_esima al tempo t, 
Vf rappresenta detta distanza, 
w rappresenta un secondo valore prefissato, e 
J indica un intomo di celle della cella i_esima; 

determinare (9) la somma in modulo delle variazioni del valore di ogni 
cella tra il tempo ^ ed il tempo /+i; 
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ripetere (10) detta fase di determinare il valore che assume ogni cella di 
detta matrice tridimensionale di celle se detta somma h maggiore di un terzo 
valore prefissato. 

2. Metodo in accordo alia rivendicazione 1 caratterizzato dal fatto che 
detta cella di detta matrice tridimensionale e di forma cubica. 

^ 3. Metodo in accordo alia rivendicazione 1 caratterizzato dai fatto ctie 
detta cella di detta matrice tridimensionale di celle e di forma parallelepipeda. 

4. Metodo in accordo alia rivendicazione 1 caratterizzato dal fatto che 
detta distanza h determinata mediante la seguente formula 

dove 

rappresenta le coordinate del centro della cella i_esima, 
p rappresenta le coordinate di im pimto, di detta pluralita di pimti di 
detto oggetto, a distanza minore dal centro della cella i_esima, 
a rappresenta un quarto valore prefissato. 

5. Metodo in accordo alia rivendicazione 4 caratterizzato dal fatto che 
detto quarto valore prefissato e compreso tra 1,5 e 2,5 e preferibilmente pari a 
2. 

6. Metodo in accordo alia rivendicazione 1 caratterizzato dal fatto che la 
fase di fissare il valore di alcune celle di detta matrice tridimensionale di celle 
ad un primo valore prefissato comprende la fase di fissare al valore +1 tutti i 
valori delle celle poste ai bordi di detta matrice tridimensionale di celle e 
fissare al valore di -1 il valore delle restanti celle. 

7. Metodo in accordo alia rivendicazione 1 caratterizzato dal fatto che 
detto secondo valore prefissato h compreso tra 0,1 e 0,9. 
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8. Metodo in accordo alia rivendicazione 1 caratterizzato dal fatto che 
rindice j rappresenta un intomo di celle a contatto di faccia o di spigolo o di 
vertice dalla cella i_esima. 

9. Metodo in accordo alia rivendicazione 1 caratterizzato dal fatto di 
comprendere la fase di filtrare detto valore che assume ogni cella di detta 
mjtrice tridimensionale di celle moltiplicandolo per xm quinto valore 
prefissato. 

10. Metodo in accordo alia rivendicazione 1 caratterizzato dal fatto di 
comprendere la fase di moltiplicare ciclicamente lungo i tre assi spaziali detta 
distanza per un coefificiente moltiplicativo. 

11. Metodo in accordo alia rivendicazione 1 caratterizzato dal fatto che 
la fase di definire le coordinate di ima plurality di punti di detto oggetto 
comprende la fase di suddividere detto oggetto in almeno due parti distinte; 
effettuare tutte le fasi rimanenti della rivendicazione 1 per ognima di dette 
almeno due parti distinte; combinare detto valore che assume ogni cella di 
detta matrice tridimensionale di ceile di dette almeno due parti distinte 
scegliendo il valore minore nel caso di sovrapposizione di celle. 

12. Programma per calcolatore comprendente un codice di programma 
che esegue tutte le fasi di qualsiasi rivendicazione precedente quando detto 
programma h eseguito su detto calcolatore. 

13. Progranmia per calcolatore registrato su un supporto utilizzabile da 
detto calcolatore per controllare Tesecuzione di tutte le fasi di qualsiasi 
rivendicazione precedente. 
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